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Abstract 



One of the greatest problems of standard cosmology is the Big Bang singu- 
larity. Previously it has been shown that non-local ghostfree higher-derivative 
modifications of Einstein gravity in the ultra-violet regime can admit non- 
^ ' singular bouncing solutions. In this paper we study in more details the dy- 

j^ ■ namical properties of the equations of motion for these theories of gravity in 

presence of positive and negative cosmological constants and radiation. We find 
stable infiationary attractor solutions in the presence of a positive cosmological 
constant which renders inflation geodesically complete, while in the presence of 
a negative cosmological constant a cyclic universe emerges. We also provide an 
algorithm for tracking the super-Hubble perturbations during the bounce and 
show that the bouncing solutions are free from any perturbative instability. 
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1 Introduction 

Primordial inflation is one of the most well motivated paradigms of the early universe, 
for a recent review, see [T]. A large e- foldings of inflation can explain the large scale 
structures of the universe and also the origin of seed perturbations for the cosmic 
microwave background radiation. In spite of the great successes of inflation, it does 
not address one of the most important aspects of the Big Bang cosmology. For 
any equation of state obeying the strong energy condition p > — p/3, regardless of 
the geometry (fiat, open, closed) of the universe, the scale factor of the universe in a 
Friedmann Robertson Walker (FRW) metric vanishes at t = 0, and the matter density 
diverges. In fact all the curvature invariants, such as R, OR,..., become singular. 
This is the reason why it is called the Big Bang singularity problem. 

Although inflation requires p < — p/3, it does not alleviate the Big Bang singular- 
ity problem, rather it pushes the singularity backwards in time. Many authors have 
pondered on this debatable issue, whether inflation is past eternal or not [2|, [3], H], 
and the conclusion is that it is not, at least in the context of Einstein gravity as 
long as the average expansion rate in the past is greater than zero, i.e. Hav > [5]. 
The fluctuations grow as the universe approaches the singularity, and the standard 
singularity theorems due to Hawking and Penrose hold, which inevitably leads to a 
collapse in FRW geometry as long as the energy density is positive [5] (see also [3]). 

There have been many attempts to resolve this cosmological singularity problem 
classically [5] , and in the context of string theory [6] , for a review see [7] , but none has 
successfully resolved the space-like singularity in the context of a flat universe. All the 
stringy toy models which have been constructed so far seem to suffer from classical 
and/or quantum instabilities and inconsistencies, such as the presence of negative 
tension branes, closed time like curves, negative-norm ghost states, etc., see [5|[9t[TU] 
for details. 

One alternative view would be to study fully non-perturbative and ghostfree higher 
derivative actions of gravity inspired from string theory [TT]: An infinite set of higher 
derivative corrections (series expansion in a', the string tension) mostly in the form 
of an exponential are often known to appear in stringy models, such as tachyonic 
actions in string field theory [12]- [19] (for a review see [20]), Bulk fields localized on 
a codimension-2 branes [2T] and various toy models of string theory such as p-adic 
strings [22] (see also [23]), zeta strings [21], and strings quantized on a random lattice 
[23 [2S] (see also [2Z|)- These corrections appear already at classical level (i.e., at 
the tree level), as well as loop level or even non-perturbatively (e.g., as expected 



in the 1/N expansion in some Yang-Mills theory [28]). Incidentally, recently it was 
shown [21] how one can consistently perform quantum loop calculations in these non- 
local models and recover, in the context of finite temperature field theory, some of 
the novel properties of string theory such as thermal duality and Hagedorn transition. 
Thus it is interesting to explore whether similar theories of gravity, albeit only at 
the level of toy models, can provide us insights into the resolution of cosmological 
singularities. 

It is worth while to point out that like most higher derivative actions, these theo- 
ries have improved Ultra-violet (UV) behavior and can render gravity asymptotically 
free, but importantly without introducing ghost states [H] , for some earlier attempts, 
see [301 [3ll [32] , which were not very successful. Note that asymptotic freedom is re- 
quired (in addition to renormalizability) in ordinary field theory for phenomenological 
reasons (scaling in processes such as deep inelastic scattering), as well as to avoid cer- 
tain problems in the nonperturbative definition of the theory. This condition usually 
takes the form of "asymptotic safety" [SB], a more general requirement for having a 
well-defined UV limit. The theories we consider are therefore natural candidates for 
such GR extensions. Building on the work done in [11], in this paper we continue 
to study cosmological background solutions and their perturbations arising in the 
higher derivative nonlocal gravity models, and in particular whether they can resolve 
the Big Crunch/Bang singularity problem of GR robustly. Previously, non-local u 
higher derivative corrections to gravity has been studied in the context of black hole 
singularities [S^, understanding ultra-violet behavior [35], cosmological applications 
(for instance, involving dark energy [36], inflation [37], and generating density pertur- 
bations [38]), and modifications to the Newtonian potentials in the metric [39]. For 
discussions on these topics in the more general context of higher derivative theories 
of gravity, see [10] Q 

In [11] non-local asymptotically safe theories of gravity were already constructed 
and it was shown that they can admit non-singular bouncing solutions. However, it 
was not clear how generic these solutions were, whether these solutions were stable 
under small perturbations, and whether these theories contain other singular solutions 
or not. The primary focus of this paper is to shed some light on these issues, and 
also see whether they can provide phenomenologically viable cosmological models. 
For instance, can one have a geodesically complete inflationary paradigm once the 
non-local modifications to gravity are included? We know that for the success of 
the inflationary paradigm, the super-Hubble fluctuations must "freeze" , is that true 
in the non-local models or do the fluctuations decay/ grow during inflation? In 
cyclic and bouncing alternatives [12]- [18], the issue of transmitting perturbations is 
of paramount importance. Can we understand the evolution of the fluctuations in 
the concrete nonsingular bouncing set-up that our model provides? 



^Some of these models are very similar to the models we consider in our paper and contain and 
infinite series of higher derivative terms, but in the literature highly non-local operators, such as 
/(l/D), have also been considered. 

■^In parallel, string-inspired non-local scalar field theories have also been studied quite extensively 
in recent times, including several interesting applications jlTj . 



In this paper specifically we will seek a two parameter family of solutions, labeled 
by A and fi, which solves the higher derivative theory of gravity. Interestingly, this 
class of solutions demands the presence of a cosmological constant, which is no longer 
free but is determined in terms of the two parameters, and radiation. One particular 
solution for the scale factor turns out to be a hyperbolic cosine bounce with a positive 
cosmological constant, where /i = — 6A^ [11]. In principle A can be both positive and 
negative in order to provide a non-singular bouncing cosmology, but for A < 0, it is 
not possible to recover the GR limit at late times. The most interesting case is when 
A > where not only we recover the right infrared behavior of gravity, but also two 
interesting classes of solutions with positive and negative cosmological constants, i.e. 
±A. 

The case of A > is very interesting from the point of view of inflation, as it 
paves the way for a geodesically complete paradigm in past and future. There exists a 
deSitter solution as an attractor for large enough times before and after the bounce. 
This is a good news for both low scale [H] and an intermediate scale models of 
inflation [5U] which are driven by the vacuum energy of the inflaton field which can 
excite the observed relativistic and non-relativistic degrees of freedom after the end 
of inflation. On the other hand the // > 0, A < case is equally interesting from the 
point of view of providing an early phase of cyclic cosmology. In this case there does 
not exist any deSitter attractor but there are also no singular solutions. Although, a 
mere cyclic phase is a rather boring scenario as it fails to reproduce the universe we 
live in, recently it has been proposed that once interactions between radiation and 
non-relativistic matter are taken into account then we can obtain a cyclic inflationary 
paradigm where the cycles become asymmetric and the universe expands bit by bit 
to give rise to an exponential inflation over many many cycles [51], see also [52]. Both 
A > 0, A < cases will be studied analytically and numerically in our paper. 

In general, taking into account a preceding contracting phase removes the infla- 
tionary infrared divergences of quantum fluctuations [53]. Here we study the evolution 
of the super-Hubble perturbations during the bounce and ensure that they do not 
feature any classical and/or quantum instabilities. Interestingly, we flnd that on large 
scales all the modes are decaying except one which remains frozen. This frozen mode 
could be useful if we wish to evolve the perturbations either prior or after the bounce 
in order to match the observed CMB perturbations. However, we will not delve into 
any details about these issues in this paper. 

We will begin our discussion with a brief recap of nonperturbative higher derivative 
action and its equation of motion in section 2. We will discuss the background 
solutions in section 3, where we also provide an analytical example of a speciflc 
class of non-local models. In section 4, we will study various cases of A > 0, A < 
and their dynamical attractor solutions for A > 0, A < 0. In section 5, we track 
the perturbations during the bounce. Some series expansions are conflned to the 
appendices: IB. II presents a general series solution for the two-parameter ansatz, and 
IB.2| IB. 3IIB. 41 consider some special cases. We conclude in section 6. 



2 Ghost-free non-local theories of gravity 



2.1 Ghost-free higher derivative actions 



In [TT], a string- inspired ghost-free and asymptotically free higher derivative gravity 
action was introduced. In this paper we investigate their cosmology in more depth 
and generality, and suggest to use this action to study how perturbations propagate 
across the bounce. Before doing this, we shall review the construction of [TT] . 

The action discussed in [TT] takes the form 

M^ f 
S = ^ d^x y/^F{R) , (2.1) 



where Mp is the reduced Planck mass and 

oo 

F(/?) = /2 + 5^-|^i?n-i?. (2.2) 

n=0 ^^* 

Here, M^ is the mass scale at which the higher derivative terms in the action become 
important. In order to ensure that the action is ghost-free, it is important that the 
series does not truncate at a finite order, i.e. the coefficients c„ must remain non- 
vanishing to arbitrarily high orders. For the purpose of illustration, let us provide a 
specific example: If we choose 

1 (-!)"+! 



6(n + l)! 
the series can be summed up and F(R) becomes 



(2.3) 



F{R) = R-^{ )R. (2.4) 

To understand the absence of ghosts, it is convenient to re-write the action in 
terms of a scalar field ip and a Lagrange multiplier field in the form 

S = i^|d*.y=5(^*fl + ,Af;<,.-|;^-[V.(*-l)--|v.^).(2.5) 

After performing conformal transformations, Qmn -^ ^gmn, and linearizing around 
$ = 1 we find 

S = -fJd^xV^9\^R+^<Pa<p + ^Y.^c^^^-^<py (2.6) 

where ^ = e'^ ^ 1 + (p. We note in passing that while is dimensionless, ip has mass 
dimension 2. The variational equations with respect to and tp yield 

ip = 300 (2.7) 



and 



2 Z^ C" , .2(n+l) ' (2- 

n=0 ^'^* 



respectively. We can now combine these equations to obtain 

= r(n)0 = 0. (2.9) 



- gn+l 

1 - 6 2^ C„ 2(n+l) 
n=0 ^^^* 



Ghosts would appear as poles of the propagator of (/). If we can sum up the infinite 
series on the right hand side of (12.91) into the form 

r(n) = e^(°) , or r(n) = (i-^\ e^(°) (2.10) 

\ m^ J 

where 7(w) is an analytic function in the entire complex plane, then it follows [11] 
that there are no ghosts. One simple example of a ghost-free action results if we 
choose 

7(w) = w, (2.11) 

which is precisely what one obtains when the coefficients Cn are chosen as in (12. 3p . In 
case one is concerned that integrating out the field ip above changes the conclusions, 
we present alternative derivations of the ghost-free condition in the appendix \^ 

One way to see that the action given by (12.11) with the coefficients chosen as in 
(12. 3p is asymptotically free is to solve [11] the linearized field equations for a point 
source and to observe that the contribution of ultraviolet modes is exponentially cut 
off. This is a rather typical scenario in these models, the UV behavior improves due 
to the exponential "regulation" . 

The asymptotic freedom of the theory leads to nonsingular cosmologies. An exact 
nonsingular solution was constructed in [TT] for a limited class of non-local actions of 
the form (12. 2p . One of the focus of this paper is to generalize the class of bouncing 
solutions, so that we can look at a much larger set of non-local actions. Also we want 
to study the stability and phenomenological viability of these models. 



2.2 Einstein- Schmidt (ES) equations 

For this purpose we first want to simplify the modified Einstein's equations for our 
action (12.21) . We will henceforth refer to them as ES equations, as much of the work 
on such generalizations was done by Schmidt [51]. The ES equations look like 

p 



where T^^, is the matter stress tensor and the generahzed Einstein tensor G^,y, called 
here the ES tensor, has the form 

G/iu = F^R^y — -Fg^y — Fo-fj^y + g^ynFo 



+ 2 5^ [9UMa''-'R)n;. - F;v;(,(n^-'^);.)] , (2.13) 

1 

where we have defined 

oo . 



M=N 



and the symmetrization convention is such that ^(/i^) = A^y + A^^. Let us obtain the 
ES tensor, G^y, for the special case when 

F = Rn'R. (2.15) 

One finds 

0% = 2R^yU^R - 2{U^R),,y + 2g^y (vT^^R - ^RD^R^ 

n 

+ - J2 [g,u{n''-"'R{D"'-'R)n;. - (□"-'"i?);(^(n'"-^i?);.)] . (2.16) 

m=l 

Since the ES tensor is linear in F, the complete ES tensor for the more general F in 
Eg. (12. 2 p is given by 

oo 

Gf,u = Gf,y + Y_, 2(n+l) '^M'^ (2-1'') 

n=0 ^'^* 

In particular it is rather useful to look at the trace of the ES tensor. Defining as 
before 

oo 

G^g^'G,y = G + Y.T7^G\ (2.18) 

n=0 ^'^* 

where 

n 

Qn ^ ggn+i^ ^ j2 [(n"-™i2).^,(n"'"^^)''' + 2(n""™i?)(n'"i?)] (2.19) 

m=l 

gives the contribution to the trace of the ES tensor from each derivative order of 

(n + 4). 

3 Finding bouncing backgrounds 

3.1 The ansatz 

Quite remarkably infinite differential equations similar to what we are interested 
in has been studied by mathematicians, see [55], and several rigorous results exist. 



especially for linear systems. The various techniques and results developed has now 
been successfully applied to analyze cosmology involving infinite differential equations 
involving scalar fields as one encounters in p-adic [SH] and string field theory [T^-|19]. 
In particular what has been realized is that mostly the non-local infinite differential 
equations can be approximately mapped to a collection of local field theories. To 
make things concrete, consider a non-local action of the form 

S = M'^ I d^x(t)T{U)(t) (3.20) 

leading to the field equation 

r(n)(/) = o. (3.21) 

Let us assume that F is an analytic function with A^ distinct zeroes of order one, at 
points rn^ with i = 1 . . . N. Then one can rigorously show that the set of classical 
trajectories for the full non-local action is simply given by the "union" of solutions 
to the following "local" field theory equations 

(n + m2)0 = O , i = l...N (3.22) 

One can thus use this correspondence to obtain solutions to the full non-local equa- 
tions by solving only the usual field theory equations. When one adds interactions 
to the theory, then such a correspondence is no longer exact, but again one can show 
that one can obtain approximate solutions to the original non-local (and now non- 
linear) equations from local field theory solutions. This procedure was successfully 
used to obtain, not only background solutions but also in computing the spectrum 
of perturbations around the background [SS]- Our aim in this section is to develop 
a similar algorithm to obtain solutions (both background and perturbations) for our 
non-local theory of gravity. 

We are going to focus on the situations where there is only a cosmological constant 
and radiation as the matter source. We will see that this situation is similar to the 
scalar field with no interaction term (I3.2ip . Thus one expects to be able to find 
exact solutions of the non-local theory by solving only "local" field equations. There 
is however one crucial difference: gravity being a gauge theory leads to additional 
constraints which are different from the "local" analogues. This explains why we 
can get qualitatively very different background solutions in the non-local theories as 
compared to the local counterparts. This will become clearer as we proceed. To begin 
with let us provide the "local ansatz" for our non-local theories. 

In order to implement the mapping between non-local and local theory, it is im- 
perative that we start with the trace equation for gravity M^G = T, where G is given 
by (I2.18P and T is the trace of the matter stress tensor: 

oo / n 

G = -R + J2 2{'n+i) 6n"+^i? + J2 [(□"""'^);M(n"^'^)''' + 2(n""'"i?)(n'"i?) 

^* \ m=l 

p 



Since radiation is conformal, it doesn't contribute to the trace equation. In close 
analogy with the scalar case we now decide to look at the subset of solutions which 
satisfies the ansatz 

nR = XR + fi, (3.24) 

where A has the dimensions of M^ and /i has the dimensions of energy density, M^. 
Unlike the scalar case we do not however claim that the subset of solutions of all 
the ansatz of the above form gives us the complete set of solutions to the full non- 
local equations. For the purpose and scope of this paper we will restrict ourselves to 
solutions of the non-local equations which also satisfy the above ansatz. 

Now, from flX^ it follows that 

n"i? = A'^i? + A"- V • (3.25) 

One can simplify the terms in the right hand side of fl3.23p . and after some straight 
forward but tedious algebra one finds 




G = A,R + A2{R'-—] +A3, (3.26) 



where 



\ \ n+l 00 / \ \ n—1 

A \ ^ \r^.„ .. C„ / A 






n=0 ^ * ^ n=l 



M? 



^-r(A)-2i^(i,l-r(A))-^).f.., (3.27) 

n=l * ^ * / 



00 /\\" ^^ / \ 



n=0 * ^ * ^ n=2 

We can now equate the coefficients of R and the square of R to zero and require the 
remaining constant to be given by the A, to obtain solutions to the trace equation. 
This implies three conditions, which can be put into the following simple form 

-l + 6coX/M^ 

^^ = ^n~2 ^' 3.30 

3A^ — jj, 

r'(A) = ^, (3.31) 

A . -^. (3.32, 

Thus given a specific higher derivative gravity action (c^'s), we can find an ansatz 
of the form (13.241) which solves the trace equation corresponding to the non-local 



action (12. ip . f l3.3ip determines the value of A, and f l3.30p then specifies the value of 
/i uniquely. However, one has to still satisfy the additional constraint (13.32 p . Due to 
our specific ansatz, the cosmological constant can no longer be a free parameter, but 
instead is given by the equation (13. 32 p . Another way of thinking about this is that 
given a specific cosmological constant, we can only find solutions of the type (I3.24p 
for a subclass of non-local actions under consideration. In section l374l we will provide 
an illustrative example of a class of non-local models admitting bouncing solutions. 



3.2 Nonsingular bounces 

Since we have "reduced" the non-local infinite differential equations to a finite fourth 
order equation (I3.24p . this is not a very difficult problem to analyze and can be done 
rigorously for instance, using power series methods. For simplicity we will focus on 
symmetric bouncing solutions, so that the Hubble constant admits an expansion like 
the following 

oo 

H = M,Y1 ^2n+i(M,t)2"+i . (3.33) 

n=0 

By substituting (I3.33P in (I3.24p one can now systematically determine the /i^'s by 
looking order by order in (M^,t)^. Actually, since (I3.24p is a third order differential 
equation in H, we expect three independent "integration constants" . By choosing a 
symmetric bouncing solution (I3.33P we have implicitly chosen /iq = ^2 = 0. Thus we 
are still left with one unconstrained parameter in hi. 

In appendix IB. II we have derived the series expansion and in appendices IB. 21 and 
IB. 41 considered some special cases. The first few /i„'s are: 

h = -^Jfi + 6hM^{X + AhM^)) , (3.34) 

^5 = ^7^ (V + 2/ii(3A2 + llfi)M^ + 132XhlMt + SSAhfM^) , (3.35) 



h 



720M6 
1 



7 



fi{3X' + 26/x) + 6/iiA(3A^ + 73/i)M,' 



90720M8 
+ UhliUlX'^ + 24:2^i)M^ + 1742AXhlM^ + 39168/i^M,*] . (3.36) 

From the structure of the series it is apparent that each coefficient /i2n+i is uniquely 
determined by the previous ones. The highest order coefficient appears always lin- 
early, guaranteeing a real solution. Thus given a specific non-local model (and the 
appropriately fine-tuned cosmological constant) we have a one-parameter solution of 
symmetric bounces, the parameter being hi, or equivalently, the value of H at the 
bounce point. An important point to note is that although H cc hi and hence can be 
as small as one wants, the time scale of the bounce is still naturally governed by the 
M^ scale. As one can see from (13. 34113. 361) all the higher order /i„'s will naturally be 
C(l),evenif /ii ^0. 



3.3 Bianchi Identities & the energy density at the bounce 

So far we have only solved the trace equation, and one may wonder whether this is 
sufficient. The answer is yes, and has to do with the fact that gravity is a constrained 
system and satisfies the Bianchi identity. For time dependent solutions 

G = -4 A ^ Goo = ^a-^ . (3.37) 

In other words, the modified Hubble equation can (at most) admit radiation which 
being conformal is traceless and therefore doesn't show up in the trace equation. Thus 
the solutions that we obtain via our "local" ansatz provides a solution to the complete 
set of ES equations (12.121) . but in so doing may inherently assume the presence of 
some amount of radiation as well. Although, there is nothing wrong with having a 
bit of radiation, we need to determine po because po niay turn out to be negative. 
This would mean having a ghost like radiation degree of freedom, precisely the kind 
of states we set out to avoid. 

In order to compute po and p one has to look at the Gqo and Gu equations: 



Goo = FoRoo + -F-Fo-nFo--J2[Fnn''R + Fn{a^~'R)]. (3.38) 

Gu = Foi?,,-yF-Fo;,, + a2nFo + y^[F„n"i?-F„(n"-ii?)] .(3.39) 

1 

where F„'s are defined via f l2.14p . For our specific action this simplifies to 

oo 

m=n 

oo 

Fo = 2j2cmO"'RioTn = 0. (3.40) 



m=0 



For the energy density and pressure we then get, using the constraints fl3.30p - fl3.32p 
from the trace equation: 



A-2cop/M,2 r / ^ , ^^x , ^ „ • .,, , 1 - 6coA/M, 



p 

t2 



Ml ^ 3A2-P r V ""; ' J ' 12A3-4Ap 



H{H\ + 2H] +6H^H-H 



* ,.2, 



A^Al) 



1 uM„ 
P ^ 3^+^. (3.42) 

consistent with having the cosmological constant and radiation. It is easy to check 
that the above two equations imply trace equation. Another consistency check is to 
compute the continuity equation, 

p + 3if(p + p)=0, (3.43) 



which now indeed holds identically. We note that the solutions do not depend on the 
details of the higher order contributions to theory once the constraints fl3.30p - fl3.32p 
are satisfied. Namely, the coefficients c^'s do not enter explicitly into the expressions, 
but their sums have to satisfy the three conditions. Only the value of cq is important, 
though only when the cosmological constant is nonzero. 

The expression for pq, the radiation energy density at the bounce point can be 
obtained as follows: Notice because of our ansatz, at the bounce point H = H = 
and H = hiM^. Thus we have: 

3M^iX-2cofi/M^)ip-12hlMt) 
^° = 12)^^Ail • ^'-^^^ 

Again we find that as long as /i 7^ 0, then even if hi — )■ 0, "naturally" we expect 
Po ~ 0{MpMl). To have a consistent bouncing solution we must make sure that 
Po>0. 



3.4 An analytic example: hyperbolic cosine bounce &; a spe- 
cific class of non-local models 

We have seen so far how we can find bouncing solutions given a non-local action. 
For the purpose of illustration, let us now look at a specific analytic example of a 
bouncing solution which obeys the ansatz fl3.24p and was already discussed in |llj : 



a{t) = ao cosh -i / —t , (3.45) 

where A > 0. This background satisfies our ansatz with the specific parameter com- 
bination 

nR = \R-Q\^ ^ /i = -6A^ . (3.46) 

The general results fl3.32p and fl3.44p imply that we need a cosmological constant, 
and in general also radiation to support this solution 

A = hul,, (3.47) 

PO = -^{\' + 2h\Mt){l + l2\c,/Ml). (3.48) 

The last equality implies we need a negative quadratic coefficient Cq < — M^/(12A). 
In our expansion for the Hubble parameter, if we choose uj = \fX and hi = 1/2, the 
general solution (13. 34113. 36|) give us 

^W^A^-iAV + lAV-JIjAV^..., (3.49) 



which is nothing but the expansion corresponding to (13.451) . We have thus shown that 
consistent hyperbohc cosine bounces with parameter A > exist in models satisfying 
the following requirements: 



r(A) = 
r'(A) = 

Co < 



5(- 

r-1 

A 

_Ml 

12A 



6coA/m; 



* ) 1 



(3.50) 
(3.51) 
(3.52) 



Clearly, A is over-determined, but there are infinite number of different sets of c„'s 
which can satisfy the above constraints and therefore admits the background evolution 
(I3.45p . For the purpose of illustration, let us provide a one-parameter family labeled 
by To of non-local actions which satisfies (13. 50113. 521) . Consider the following r(n): 



r(n) 



exp 
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(3.53) 



By solving (I3.30p and (13.3 1|) for the hyperbolic cosine bounce, and imposing 



To = exp ai 

we can parameterize the a's as 

ai = 

^2 = 
04 = 

where we have used/chosen 
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The bounce radiation density is given by 



(3.54) 

(3.55) 
(3.56) 
(3.57) 

(3.58) 
(3.59) 



Thus, as long as Fq > 1 the radiation density is positive and we obtain a consistent 
hyperbolic cosine bounce in the nonlocal model. Below we show the values of a's 
as a function of the Fq parameter for which the solution exists. After some series 
manipulations, we find that the coefficients c„ in the action corresponding to this 
theory are given by 
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Figure 1: Values of the three parameters a as a function of Fq as indicated in the 
figure. 

where only the terms where the factorial is a positive integer are counted in the sums. 
Thus we have shown how one may also reverse-engineer nonlocal theories to generate 
the desired background evolution. 

In figure ?? we provide the range of the relevant parameters as a function of Fq 
for A = M2. 



4 Cosmologically viable non-local models 



To summarize, we have found that given a specific higher derivative theory of the 
form (12. 2p . i.e. given a specific F(n) function, we can find a one parameter family, 
labeled by the first derivative of the Hubble parameter at the bounce, of symmetric 
bouncing solutions. Each coefficient /i2n+i is uniquely determined by the previous 
ones, and thus one may extend the solution in principle arbitrarily late in time. 
However there are several important questions one needs to answer before we can 
start taking them as serious candidates for consistent and phenomenologically viable 
non-singular theories. In this section we want to discuss these issues. 

The first question has to do with the radiation density at the bounce which is 
fixed in terms of the model parameters and initial conditions (this could be turned 
the other way around too), and is given by (13.441) . Since we want to consider the- 
ories without ghosts, this radiation energy density must be positive which provides 
additional constraints on the parameters. The second issue has to do with whether 
suitable GR limits of the solutions exist? Finally, it is important to find out whether 
the bouncing solutions that we have found are generic, or are there background so- 
lutions which end up in a singularity? In this section we will try to address these 
questions staying within the confines of our ansatz equation. 



4.1 A < Case, problems recovering a GR limit 

Before we discuss the non-local theories let us revisit the much studied R^ modification 
of gravity: 

F{R) = R+^R\ (4.61) 

In Ref . [57] it was shown that a zero- measure set of nonsingular solutions exists for this 
model. Even in the absence of radiation and a cosmological constant, the cosh-bounce 
is a solution in a closed universe [5H]. Nonsingular solutions and attractor solutions 
for conformally invariant fields have been found also taken into account quadratic 
corrections to the gravity action [60], see also [61]. The bounce solution that we have 
for R^ gravity is a reiteration of these results. 

This quartic higher derivative action is actually unique because it is the only case 
when the solutions are not over determined: 

, M2 ^ 2AM2 ^^. 

A = ^ and u = *- . (4.62) 

This was checked in the appendix IB. 31 One can also check that there is a large 
parameter space where the radiation density is positive. For instance, for the special 
case when the cosmological constant and fi vanishes, the constraint (13.441) gives us 
po = —3Mphluj^/X. Thus as long as A < we can have a consistent bouncing solution. 
This is rather puzzling because it is also known that in a fiat universe, if we allow 
only a positive cosmological constant and radiation density, there are no viable non- 
singular solutions (even for slightly more general models with an additional term due 
to conformal anomaly) [59]. So, does this mean that just the simple R^ modification 
can avoid the Big Bang singularity problem after all? 

The answer turns out to be rather interesting. Let us try to understand the 
attractor behavior of the ansatz equation when A < 0. Let us start with a late time 
trial solution 

Hit) = hit. (4.63) 

(EMD then reads 

-12h\t^ = I2\hlt^ + Oil) 

Thus we have a late time super-infiationary attractor solution given by 

hi^-^. (4.64) 

6 

In the expanding branch, since /ii > 0, this means that the attractor solution only 
exists if A < 0. Our numerical studies confirm this behavior, see figure [3], and we will 
also provide a "Dynamical System Analysis" further corroborating it. 

The upshot is that as long as A < 0, although these models admit nonsingular 
bouncing solutions, they are phenomenologically not viable because we cannot recover 



the GR phase but rather are stuck in the QG phase foreveio. Hence, from now on 
we are going to focus on non-local theories with positive A. The difference to pure 
R^ gravity is clear from Eq.f l3.24]) . This is the trace equation for i?^-gravity, but 
in that case both A and fi are uniquely fixed by fl4.62p . and one inevitably obtains 
negative radiation energy density for the nonsingular solutions. In a nonlocal theory 
of the form f l2.ip . cq of course doesn't alone fix these coefficients, and hence we can 
get bouncing solutions with just positive-density sources. 



4.2 A > case, and the issue of singularities 

We have seen that if A > the QG super- inflationary late time attractor no longer 
exists paving the path to recover GR cosmology at late times. However, we still have 
to worry about singularities? Having a single bounce does not guarantee that the 
evolution encounters no singularities in the future (or the past). Can we find if and 
when we have completely nonsingular evolutions? In order to answer this we will look 
at the two cases of positive and negative cosmological constant separately. 



4.2.1 A > 0, Geodesically complete inflation 

For infiationary applications we will be interested in the case when A > 0. The 
expression for the cosmological constant f l3.32p then tells us that /i < 0. From fl3.44p 
it is now clear that to have a positive po 

A - 2co/i/M2 < . 

Since p < and A > 0, this means Cq < 0. Also note that we can always choose 
I Co I = 1 by choosing M* appropriately. Hence setting cq = —1, the inequality is given 

by 

.<-^. (4.65) 

Now, since our universe is filled with radiation and vacuum energy, for a positive 
cosmological constant the GR solution corresponds to a deSitter solution. It is obvious 
that one of the late time solutions is the deSitter universe: 

The important question we want to find out is how generic and stable this late time 
deSitter solution is. Both numerical and analytical studies confirm that this is indeed 
a stable attractor, we will see this in the next subsection 14.31 However, are there 



^However, when considered within the first order formalism, which is an alternative way to avoid 
some problems in higher (finite) derivative extensions to gravity, the quadratic model of the form 
(|4.6ip does generate a bounce with a smooth transition to the GR phase [5^ . 



other singular attractor solutions which may sometimes prevent the universe to reach 
the above deSitter background? 

Consider the ansatz 

H{t) = ^ (4.67) 

Then as t — t- 0, the most singular contribution in the ansatz equation (I3.24p is (9(l/t^). 
If we demand that the coefficient of the 0{l/t^) term must vanish we end with two 
possibilities: 

p = - or p = 1 . (4.68) 

The ffist one is the usual GR solution in radiation dominated era. This is in fact an 
exact solution of (I3.24p when A = /i = 0, but crucially this solution is unstable, as 
we will discuss in the next subsection. The second solution however represents a Big 
crunch singularity. Can these singular solutions be projected out when we restrict 
ourselves to having positive radiation energy densities? A simple way to check this is 
to compute the leading singular behavior in pr as we approach the singularity. From 
fl3:4T]) we find 

Pr \ + 2ii/Ml ^^/l 



(I4.65P then implies that radiation is indeed positive. Thus we cannot exclude the ex- 
istence of singular solutions and they represent genuine attractors. However, as will 
be shown in detail below, for symmetric bounces we are considering, the singularity is 
never reached, and the universe evolves asymptotically towards the global GR attrac- 
tor solution at low curvature. Highly asymmetric bounces however can sometimes 
end in a singular Big Crunch singularity, but there is a large region in the parameter 
space of initial conditions for which we reach the deSitter phase asymptotically at late 
times. This will be discussed in more details in section 14.31 Finally, one may worry 
about other higher order singularities {H oc l/t" with n > 1), but simple counting 
arguments show that such singularities don't exist. 

Incidentally, using (I4.65P one can also obtain a lower bound on the value of the 
cosmological constant. We have 

A = -^^^>_^^. (4.70) 

4A - 8 ^ ' 

Note that we can adjust the cosmological constant to lower values by lowering the 
scale M* (or equivalently, by keeping M* fixed but allowing higher quadratic term cq 
if one doesn't choose the convention |co| = 1 we used in this subsection). 



4.2.2 A < 0, and cyclic universes 

A rather interesting case emerges if one considers A, // > 0. In this case there are no 
QG super-inflating attractors, and neither do the singular solutions exist for positive 
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Figure 2: Cyclic evolution for some parameter combinations a,/i > 0. The Hubble 
rate H and the time t are in units of M^,. 



radiation energy densities, see Eq.f l4.69p . What we obtain instead are cyclic universes, 
where QG effects trigger a bounce, while the universe turns around once the radiation 
energy density dilutes and the negative cosmological constant cancels it. In figure [2] 
we show numerical solutions of such an evolution. 

What is especially nice is that there is a special class of non-singular actions with 
A, yU > 0, and fi < 3A^ which does not admit any singular solution. This is because 
f l4.69p tells us that singular solutions can only exist for negative radiation densities. 
Conversely, if we have standard radiation, such solutions do not exist. 

Phenomenologically, if we just have a negative cosmological constant, it is not 
interesting because we end up in an eternal cyclic scenario where the period is too 
short for structure formation. However, recently similar scenarios have been discussed 
where once one takes into account of interaction between different species and entropy 
production, viable cosmological scenarios emerge [51] . 



4.3 Dynamical system analysis 



The asymptotic state of the system can be studied by identifying the attractors in 
the phase space. In terms of the dimensionless variables 






PR 



3M^H^ 



S = 



where we have defined the auxiliary quantity 

A-2co/i/M2 



e = 9- 



3A2-/i 



^H^ 



(4.71) 



(4.72) 



the ansatz, Friedinann and continuity equations can be written as 

S' = 2eS, (4.73) 

n'j, = 2{e-2)nn, (4.74) 

e' = ^(e-2)e+^(eA-3fi^ + ^eVS)S, (4.75) 

e" = [7e' + 6(2 - e)e] (e - 1) + (2 - e)eAE + ^eVs' . (4.76) 

The system has the following fixed points: 

• A: The radiation dominated solution e = 2, S = 0. The linearized matrix for 
perturbations about this fixed point has two positive eigenvalues, meaning the 
solution is unstable. 

• B: An infiating attractor S = e = Qji = 0, which is stable. 

• C: The deSitter solution S = — 12A/(^/i), e = Qr = 0. It is an attractor when 
// < 0, otherwise a saddle point. 

To understand the meaning of these fixed points, it is useful to consider the 
expansion-normalized energy contributions 

The first two are the relative contributions from radiation and the cosmological con- 
stant, respectively, and the latter two represent the non-local corrections to gravitjQ. 
In terms of these, the Friedmann equation is 

1 = sgn{^X)QR + ^a + Qi + Qs • (4.78) 

The sign appears in front of ^2^ since its contribution can be negative if .^A < 0, 
though we always consider positive physical radiation energy density. The unstable 
fixed point A corresponds to Qf> = 1, the attractor B to Qi = 1 and the fixed 
point C to Q\ = 1. It is then clear that B is the solution (14.63114.641) we derived 
previously. Since we are interested in the post-bounce evolution, we consider the 
expanding branch, where the fixed point B only exists for negative A. Then, the 
universe will be super-infiating, but approaching e — ?■ instead of ending in a Big 
Rip. This possibility would be interesting in the dark energy context, but in the 
present study we do not regard evolutions ending in B as viable trajectories, since 
the phantom expansion is driven by Qi and we want to recover pure GR. 

Thus, to underpin our reasoning in the preceding subsections, we have shown that 
we need a positive A to avoid the attractor B and a negative /i to be attracted towards 



^Despite this intuitive form, these turn out to be less suitable for dynamical system analysis, 
since the phase space spanned by these variables contains only the fixed point C. 
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Figure 3: Dynamics at the bounce for two models for normalized relative energy 
contributions, Qi being the nonlocal and Qr the radiation contribution. Left panel: 
evolution in the A < 0, /i < case: the super-inflating solution B is reached after the 
bounce. Right panel: the same model but with flipped sign of A. Now the de Sitter 
C attractor will be reached, and both Qi and Qr vanish asymptotically. 

the GR evolution C where the nonlocal terms Qi and Q2 both vanish. Example 
evolutions of bounces evolving to B and to C are shown in Fig. [3l Numerical 
experiments confirm that the GR phase C will be generically approached after the 
turnover, given just the correct sign choices. The cosmological constant need not to 
be large for this, but can be adjusted to the desired scale. Evolutions for both large 
and small A-terms are shown in Fig. HJ Moreover, the convergence to the attractor 
is insensitive to small perturbations in the initial conditions. Since H vanishes at 
the bounce and we can normalize H to unity in units of M*, the most general initial 
conditions are given by varying H at the bounce, or equivalently, considering nonzero 
/i2- Only if /i2 is large and negative, we fail to recover the GR regime, as illustrated 
in Fig. |5]for several choices of A. 



We remind that we have only exhausted the possibilities within ansatz we have 
made. In reality the dynamics will change as matter becomes non-relativistic and 
dust sources render our system based on the ansatz invalid. Furthermore, even when 
the energy sources consist only of radiation and A, there can be solutions which 
do not obey equations ( 13. 30113. 32]) . i.e. not each of the contributions to the trace 
equation vanishes identically. Our stability analysis here concerns only the reduced 
phase space where the trajectories obey these conditions and the ansatz fl3.24p . In 
principle this class of trajectories need not be stable among the space of all solutions. 
At a more fundamental level, ghost-free condition of section (12. ip is only about the 
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Figure 4: Convergence to the deSitter attractor C when A = 1. In the left panel 
/i = —0.1, and in the right panel fi = — e~^°. The oscillations take longer for smaller 
/i, but the GR attractor C is reached eventually. 



-0.5 



-1 



h 



2-1.5 



-2 



-2.5 








^^•""'^'--^^..^^^^ 




\ t » • 




\ * * 






x=^o 


---x=^ 


I 1 


■--?i=0.1 




?i=0.01 


> 

1 
1 I 



10" 



A 



10^ 



Figure 5: Lower limit on the initial perturbation /12 to converge to the GR attractor C 
when units are such that hi = 1 ai the bounce. Above the lines, we reach the deSitter 
fixed point C. Below the lines, with large negative H (presumably an unnatural 
boundary condition at the bounce), we are driven to the big crunch. This shows that 
the GR attractor C is stable to small perturbations. 



Minkowski space, which may not be enough to guarantee that ghosts wouldn't appear 
at a nonhnear level in strongly curved backgrounds. In a complete analysis, one would 
have to include perturbations in the metric as well and allow inhomogeneous modes. 
To fully address these problems, one has to device other techniques to solve the full 
non-local ES differential equations. This is beyond the scope of our present paper. 



5 Perturbations 

In this section, we show how the non-singular bouncing solutions satisfying our ansatz 
can be used to study the evolution of large-scale perturbations by either series methods 
or by numerically solving differential equations. 

5.1 Tracking perturbations across bounce 

We want to investigate how we can track the evolution of small perturbations around 
a specific bouncing background solution. This is obviously a hard problem and we fall 
back on our intuition that was developed by dealing with scalar fields. In the scalar 
field case perturbations around a background satisfied the same local differential 
equation as the background. Since we are looking at small perturbations around the 
background solution which satisfies a "local field equation", it seems likely that the 
perturbed solutions also satisfy the same "local field equation" , by virtue of which 
they will also solve the full non-local equations. This is where the covariant nature of 
our ansatz comes into play. This strategy was successful in the scalar field case and for 
our gravity theory, in analogy, all we need to do is to perturb our ansatz fl3.24p . In this 
approach, at worst, we might be throwing away some modes of perturbations (which 
solve the full ES equations f l2.12p . but not the ansatz equation fl3.24p ). but at least 
we should be able to successfully track the ones which satisfy the ansatz equation. A 
further drastic simplification arises if we assume that there are no more than ~ 70 
e-foldings of inflation after the bounce. In this case, the physical wavelengths of the 
perturbations, that we see today in the CMB sky, at the bounce point are still too 
large as compared to the Planck length. In other words, the perturbations behave 
as super-Hubble modes during the bounce, the time scale of bounce is too short as 
compared to the physical wavelength. We can therefore ignore any spatial dependence 
of the perturbations and just look at the fluctuations of the zero- mode or the scale 
factor. (Note, if the fluctuations in the metric just depend on time, one can always 
perform a coordinate transformation such that the only independent function is the 
scale factor.). Thus, just by looking at time dependent fluctuations of the ansatz 
equation we can unambiguously track the perturbations across the bounce, and this 
formalism will be explained in details in [63] . 

Now fluctuations of the scale factor, 6{t), can be encoded as follows: 

ds"^ = -dt^ + a^{t)dx^ with a{t) = d{t)[l + 5{t)] , (5.79) 



so that 

H{t)=H{t) + 5. (5.80) 

The hnearized expressions for R and Di? is given by 

R = R + 6{6 + AH6)=R + 6R, 
nR = nR-{d"R + 3H6'R + 36R). (5.81) 

Thus the perturbation equation reads 

6"R + 3H6'R + 36R=-X6R. (5.82) 

This is a hnear fourth order differential equation and in principle can be solved. In 
practice, one has to again employ power-series methods to solve for 6(t). However, it 
is clear that a trivial solution of the above equation is 5 = 0, or 5 is a constant. This 
is analogous to the statement in GR that the super-Hubble perturbations remain 
constant. However, a crucial thing one has to check is whether any of the above 
modes are spurious or not. In other words, one has to substitute 6{t) obtained by 
solving (I5.82P in the Gqo equation and check that it is also satisfied, just as we did 
for the background. 

Let us then expand S as 

6{t) = 6i{Mj) + 62{MjY + 63{Mjf + ■■■ (5.83) 

From zeroth to the third order, the coefficients in the expansion of the perturbed 
ansatz equation (15.821) as a power series in t then give 

= 62\ + i{262hi)M^, (5.84) 

= 76ii2 + 18{6ihi - 26s)XM^ - 12(60^5 + hi{4563 - 26ihi))M^ , (5.85) 

= llS2fi + 6i7S2hi-6Si)XM^-24:[A5Sfi + hi{33di-252hi)]w^, (5.86) 

= --6i\fi + i [I5663fi - 5i/ii(9A2 + 29fi)] Ml - 6 [20^5 + hi{75ihi - 4053)] AM,^ 

- 8[63057 + /ii(435(55 + 8/ii((5i/ii-653))] . (5.87) 

where we have used the solutions (13. 34113. 361) for the background quantities in terms 
of hi. From the structure of the expansion we see that once 82 is given, each 52n 
where n > 1 is determined uniquely. Since 52n appears linearly in the coefficient of 
the t'^i"'-'^) term, the solutions are guaranteed to be real. Similarly, given 5i and ^3, 
each 62n+i where n > is determined by the constraint from the t^n-s coefficient. To 
fix 62, we can use the Friedmann constraint at the zeroth order. This gives 

3(A - 2cofi/Ml) [5ofi - 12h{52 + 6ohi)M^] = . (5.88) 

This allows us to solve 62 in terms of 6q (barring the special case M^X = 2co/x, 
when both 60 and 62 will be left as free parameters). It follows that all 62n will be 



proportional to 6q, which may be regarded as an initial condition. Up to the seventh 
order, the full solution is then explicitly 

S2 = So (-hi + ,^ J , (5.89) 

\ l2hiMtJ ^ ' 

^' = ^' luh^m ' ^ ^^ 

75ifi + ISiSih - 26s)XM^ + 12hi{2Sih - ^5Ss)Mt . ^^ . 

^' = 720Mt ' ^'-^'^ 

^' = ^° 12960/.,M8 ' (^-^^^ 



S7 



106iXfi + 3 [12(^3 - 6ihi)X^ + (104^3 - 87Sihi)fi] M^ 



30240M6 
+ 42hi{7263-196ihi)XM^ + 3Qhi{49963-306ihi)M^ . (5.93) 



This could be continued to track the evolution further from the bounce at t = 0. Only 
the constant part of the Friedmann constraint, 05. 881) contained nontrivial informa- 
tion. The solution 05.89115.9^ satisfies the time-dependent part of the Friedmann 
constraint identically. We also note that as for the background, the details of the 
higher order derivatives, i.e. c„ where n > 0, are not important for the perturbations. 
This is because the superhorizon modes have the same symmetries as the background. 
Inhomogeneous perturbations could in principle be sensitive to the higher order terms 
in the action. 

To summarize, once the background evolution is known, by fixing three initial 
conditions the perturbations can be determined order by order by the method above. 
Alternatively, one can integrate the differential equation 

(5(^) + 7(H + d) 6'^^^ +4('6 + 6H6 + 36^ + 2H + 2H^ + \\5 + 2{6H + 66 + X)S^ 
+ (7H + 4H{6H + X))6 = 0. (5.94) 



where we have omitted the bars from background quantities. In the linear approxi- 
mation this simplifies to 

5(^) + 7H6^^^ + {8H + 12H^ + X)6 + {jH + 4H{6H + X)) 6 = . (5.95) 

Below we provide analytical and numerical solutions to these equations. 

5.2 Example of hyperbolic cosine solution 

A more complete analysis of perturbations will be presented elsewhere, but here we 
look at a special case to understand how the perturbations evolve. Let us go back to 
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Figure 6: Numerical solutions of the perturbation evolution in the cosh-bounce for 
different initial conditions at tj = 0. Units are such that M* = 1. Solid (black) 
line: 6i = 6i = 1.0e~^. Dashed (red) lines: 6i = 6i = 6i = 6\ = 1.0e~^. We also 
plot a decaying mode with dash-dotted (blue) line: 6i = 6i = — 1.0e~^. Note that 
the linear perturbation 6 can always be rescaled as the normalization does not affect 
the evolution. The results confirm that the superhorizon perturbations freeze in the 
expanding phase. 



the hyperbolic cosine bounce, so that (I5.95P reads 



S^'^ + 7J^tanh(J^ti5(^) + 



7A^ 
A 



tanh 



-f- 4\/-tanh 




t sech 




6A tanh 




-t j I SAsech^ I \ -t 1 + A 



6 = 0. 



(5.96) 



The above equation can be solved numerically. We also show numerical examples of 
the perturbation evolution in figure [6] for three choices of initial conditions. As the 
perturbations are small, the solutions of both the linear equation (15.95^ and nonlinear 
equation (15.941) agree (we have checked that in the case of large perturbations, the 
effect of nonlinearities is to stabilize the growth). We see that there exists a linear 
growing mode in the contracting phase. Consequently, the amplitude of the pertur- 
bations is enhanced before they freeze in the expanding phase to a constant value, 
precisely what one would like for the inflationary mechanism to work. 



One can actually understand this feature analytically as well. At late times, the 



equation simplifies considerably: 



5(^) + 7W|5(=^) + 7X6 + ( 4 J^ ]5 = 0. (5.97) 



The solutions are of the form 6 = exp(z/t) with 




z/^ + 7W -z/^ + 7Xu^ + [ A\l '— ] u = 0. (5.9J 



{0,-1/2,-2,-1} (5.99) 



The solutions are 

V2A 

Thus all the solutions are decaying except the first one which is the usual constant 
mode, and the standard inflationary mechanism can operate successfully. 



6 Conclusions 

In this work we have shown how ghost and asymptotically free non-perturbative ef- 
fective actions of gravity motivated from string theory are able to realize non-singular 
bouncing cosmologies. We studied a class of solutions with two parameters, which 
solves the non-local and non-linear equations of motion for gravity in presence of a 
cosmological constant and radiation energy density. In particular, a specific class of 
non-singular bouncing solutions is given by a hyperbolic cosine scale factor of the 
universe. 

We studied cosmological scenarios with both positive and negative cosmologi- 
cal constants. In the former case we have found that an inflationary (deflationary) 
solution arises as a late time attractor after (before) the bounce, the non-local contri- 
butions to gravity vanishes, and the space-time asymptotes towards the GR deSitter 
limit. This result is very relevant for inflationary cosmology as now time-like and null- 
like trajectories can be made geodesically complete within the inflationary framework. 
Without the non-singular bounce the trajectories in the past light cone would always 
be flnite, which would render inflation geodesically past incomplete. Furthermore, 
the Bunch-Davies vacuum would be pathological in the infrared. 

When the cosmological constant becomes negative an interesting stable solution 
arises which leads to a cyclic universe in presence of radiation. Although just with 
radiation a cyclic universe does not lead to any phenomenologically viable cosmology, 
exciting cosmology can arise if we add a non-relativistic matter component and inter- 
actions between radiation and matter. As shown in [^I], such an interaction would 
lead to a cyclic inflationary scenario with distinct signatures for the cosmic microwave 
background perturbations. Also, our non-local model provides a concrete setting to 
study how perturbations are transmitted through a nonsingular bounce which is im- 
portant for cyclic and bouncing models in general. Finally, the non-local actions with 



A > 0, A < and fi < 3A^, constitutes an especially nice class of actions because 
they don't seem to contain any singular solutions (at least within the confines of our 
ansatz), which is not the case for actions with A,A > 0. An important task for the 
future is to determine whether these non-local actions admit any singular solutions 
at all or not. 

In order to make sure that both infiationary and cyclic attractors are stable, 
and there exists no instabilities during the bounce, we studied the long wavelength, 
super-Hubble perturbations. We found that all the modes of the long wavelength 
perturbations are decaying except the one which is frozen. The presence of a constant 
mode is a good news for seeding density perturbations for the cosmic microwave 
background radiation during, before or after the bounce. The perturbation analysis 
suggests that there are no pathologies at least on large scales which would affect 
the non-singular bounce solution. Studying the sub-Hubble perturbations is more 
difficult in higher derivative theories of gravity. However, this is necessary to consider 
the ghost problem in the bouncing backgrounds, whereas we have yet only shown that 
the theories are not plagued by this problem in Minkowski backgrounds. We hope to 
revisit this issue in a separate publication. 

To summarize, we have found non-singular bouncing solutions for the scale factor 
of the universe in a special class of nonlocal gravitational models which are ghost and 
asymptotically free in the ultra violet, but nevertheless recover the GR limit at large 
distance and time scales. We also made preliminary progress to study perturbations 
around the bouncing background and in particular demonstrated that the bouncing 
solutions are late time attractors and perturbatively stable. 
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A Condition for "ghost-freeness" 

In section 12. ![ we explained how the condition for ghost-freeness arises by looking at 
the field equations and identifying the poles in Eq. fl2.9p . An equivalent, and perhaps 
a more rigorous, way to see this is to work in terms of the action involving the two 
auxiliary scalar fields and ip directly. On a Minkowskian background the momentum 
space action looks like 

S = -^ I d^k<i/klC{-k^)^k, (A.lOO) 



where we have defined 



m3 



Now, as long as the determinant of the kinetic operator has at most a single root 
we can be assured that there are no ghosts in the theory. One can see this in the 
following way. One can perform a linear transformation to diagonalize the kinetic 
operator so that it would be of the form 



/C' - 1 „ ^ , ,^2\ 



Ti{-e) 

r2(-F 

It is now clear that the number of perturbative physical degrees of freedom simply 
corresponds to the number of zeroes of the product ri{—k'^)r2{—k^) which is nothing 
but the determinant of the kinetic operator. This in turn means that we can simply 
count the zeroes of the original kinetic operator /C(— /c^) to determine the perturbative 
states. As long as we have a single or no zeroes, we are perfectly safe from any 
problems with ghosts. For more that two zeroes there is always a ghost while the 
case with two zeroes may or may not have a problem depending on the sign of the 
residues at the poles. We do not explore this possibility here. 

Now, it is easy to check that 

det(/C) = -^r{-k^) (A.102) 

and thus indeed we obtain the same condition for the action to ghost free as deduced 
from 12.91 It is also straightforward to show that the same condition emerges if one 
integrates out the field and considers the propagator for the field ip. We then obtain 
the equation of motion 

n-^T{n)ip = o. (A.103) 

In the special case /c = it is not legitimate to integrate out 0, and when 
S^o '^ri ( -j^ ) = one cannot integrate out ip. However, as we saw above also 
in these special cases the same condition for F determines the nature of the propa- 
gating degrees of freedom. 



B Series expansions 

For the matter content we assumed a A-term and conformal matter, and for the 
metric the ansatz f l3.24p 

nR = \R + n. (B.104) 



Then the Friedinann equation fl3.4ip was derived, 



M2 
p 



9 



A - 2co/i/M2 



3A2-/i 



6coA/M2 



12A3 - 4A/i 



;^.105) 



H (HX + 2H] + eiJ^if - H^ 
To look for bouncing solutions, the series expansion (I3.33p for H was made: 

oo 

H = M,Y1 /i2n+i(M,t)2-+i . (B.106) 

n=0 

It follows that the scale factor behaves as 

CL[0) Z o 4o 

{h'l + 12hlh + 12/^3 + 32/ii/i5 + 48/i7)(M,t)^ + C((M,t)^°) (B.107) 



384 



and in the following we set a(0) = 1. The energy density is then expanded as 



Pit) 



' 4A 



+ Po 
1 
+ 6 



1 - 2ht(M,tf + (2ft? - fts) (Af.tf - 5 (2ft? - Sftifts + fts) (MJ)' 
{4ht - Uhfh + 3/^3 + 8hih - 3/17) (MJ)^ + 0((M,t)i°)j .(B.IOS; 



The first term is the cosmological constant given by the constraint (13.301) . and the 
rest of the terms describe the radiation density. 



B.l Solutions for the Hubble parameter 

In the following we assume 3A^ 7^ fi. To solve for the coefficients /i2n+i we use the 
expansions for the matter ( JB.IOSP and for the Hubble parameter (JB.106P in the two 
equations f lB.104p and f IB.lOSp . Order by order, the ansatz equation ( IB.104P yields 
the following constraints: 

= fx + 6hi\M^ + l2{2hl + 36h3)M^, (B.109) 

= {2hj + 3h3)\ + 6{2hl + llhih3 + l0h5)M^, (B.llO) 

= (4/ii/i3 + 5/i5)A + 6(10/i?/i3 + 13/^3 + 30/ii/i5 + 35/i7)M2 . (B.lll) 

Leaving M* and hi as free parameters, as they correspond to the initial conditions 
for the Hubble parameter and its first derivative, we can solve the three following 
coefficients as 

h = -,^(;i + 6/iiM2(A + 4/iiM,2)) , (B.112) 



h, 



36M2 
1 
720M6 
1 



(A^ + 2/ii(3A2 + llfi)M^ + l32XhlM^ + 384/i?M,^) , (B.113) 

/i(3A2 + 26/i) + 6/iiA(3A2 + 73fx)M^ 

(B.114) 



90720M8 
12hl{UlX^ + 2A2^i)M^ + 17A2AXhlM^ + SdlGShJM: 



From the structure of the series it is apparent that each coefficient /i2n+i is uniquely 
determined as by the previous ones. The highest order coefficient appears always 
linearly, guaranteeing a real solution. Let us then consider the Friedmann equation 
(IB.IOSP to determine the density po- The ffist three nontrivial constraints are 

= 3A/i-6co/iVM,'-4(3AWA^p+/^)Po/M2-36/i?(A-2co/i/M2)A#,115) 

= (3A2-p)2/iipoA'i'p+9/ii(A-2co/i/M2) [;/iiA + 6(/i? + /i3)M2]AO,116) 

= -po(2/i? + /i3)(3A2-p)/Mp2 

+ 9(A - 2cofi/M^)M^ [3 {lOhhj + lO/ig/ii + hi) Ml + 2hih-i\] . (B.117) 

By construction, we expect all of these to be identically satisfied, once the solution to 
the ansatz equation is used in conjunction with one of them to fix po- This is because 
the form of Friedmann equation was already dictated by the Bianchi identity given 
the ansatz we have. Thus the redundant system flB.115tlB.117p serves as a consistency 
check. The density we obtain is 

3Ml{\-2c^^/Ml){^^-l2hlMt) 

p' = m^^i;^ • ^^•''') 

The parameters of the model should be chosen such that p > 0. 



B.2 Special case fi = 3A^ 

The hyperbolic cosine bounce a{t) = qq cosh w |t belongs to the class of models 

satisfying our ansatz, in particular when p = — 6A^. There are solutions of this type 
and its generalizations admitted by the nonlocal field equations as we saw above. One 
special case which our treatment in ( IB. II) however does not cover is when p = 3A^. 
Let us consider now this. Now the constraints f l3.30tl3732l) dictate we should have 

^ = i, A = -2aM,^ (B,119) 

So the quadratic coefficient is fixed by A, but now the propagator r(A) is not deter- 
mined by the trace equation and is left as free variable. We cannot now write the 
energy density as Eq.( lB.105|) but obtain instead 



^ = 3(i-r(A)) 

p 



H^ + \ (gH^H + 2HH - iJ2 
A 



^r(A)A. (B.120) 



The three first conditions from the ansatz equation (IB. 1041) are equivalent to 
Eqs. (lB.109IIB.lll]) . The ffist one can be now simplified to 

= A^ + 2hiXMl + 4{2hl + 3h)Mt , (B.121) 



the rest remain the same. From the constant part of the Eq. flB.120p we get the 
expression for radiation density as 

?,Mi{i-nm\^-Ah\Mt) 

Po = ^ ■ (B.iiij 

Again we can check that the constraints from the Friedmann equation expansion at 
higher order in {M^,t) are identically satisfied. 



B.3 Quadratic model 

To briefly check the consistency of the solution for the model (14.611) that includes 
only the near-trivial truncation of the theory, we recall that in fourth order gravity 
the Friedmann equations can be written as 

3f{H + H^)-h = 3Hf-p, (B.123) 

f{H + 3H^)-h = f + 2Hf + p, (B.124) 

where / = dF/dR. If we now plug in f l4.6ip and say H(t) = hiuH + . . . , at ze- 
roth order the above equations yield p = lSM'^CQh\M^ / a^ and p = 6MpCo/i^M^/a^, 
matching with our constraints. The expansion of H, which does not terminate, could 
be continued to higher orders. Thus solutions can be obtained, but they are not 
cosmologically viable as shown in 14.11 



B.4 Harmonic Ricci-spaces 

Also the case of vanishing A is a special case apart from the previous ones considered in 
IB.l|B.2t since we have used the assumption A 7^ there. Now the constraints (13.301 
I3.32P from the trace turn out to imply the following relations between the model 
parameters 

ci 1 . 3/iCo fi^C2 mioc;^ 

Mj = 2^' ^ = -2Mf-2M!' '' = '■ (^-^25) 

These are impossible to satisfy. So in our case we must have nonzero A. However, a 
different set-up can be considered when the trace equation is T = —Ra.s usually and 
not identically satisfied. The simplest example is /i = A = 0. In these spacetimes the 
solutions of GR and the nonlocal model coincide exactly. Perturbations can of course 
behave completely differently, and therefore the stability of these spaces would be 
interesting to study. However, since these spaces are not bouncing we don't pursue 
this further here. 
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